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We combine experiments with simulations to investigate the fluid-structure interaction of a flexible 
helical rod rotating in a viscous fluid, under low Reynolds number conditions. Our analysis takes 
into account the coupling between the geometrically nonlinear behavior of the elastic rod with a 
non-local hydrodynamic model for the fluid loading. We quantify the resulting propulsive force, as 
well as the buckling instability of the originally helical filament that occurs above a critical rotation 
velocity. A scaling analysis is performed to rationalize the onset of this instability. A universal 
phase diagram is constructed to map out the region of successful propulsion and the corresponding 
boundary of stability are established. Comparing our results with data for flagellated bacteria 
suggests that this instability may be exploited in nature for physiological purposes. 


Bacteria often rely on the deformation of filamentary 
helical structures, called flagella, for locomotion [T]. The 
propulsion arises from a complex fluid-structure inter¬ 
action (FSI), between the structural flexibility of the 
flagellum and the viscous forces generated by the flow. 
This FSI may lead to geometrically nonlinear deforma¬ 
tions mm, which in turn can be exploited for turn¬ 
ing a, tumbling [S], bundle formation and polymor¬ 
phic transformations izli- 

Resistive force theories (RFT) [9l[T0] are often used to 
model the role of viscous forces on flexible filaments m, 
at low Reynolds number. These simplify the viscous 
loading by introducing local geometry-dependent drag 
coefficients. More sophisticated descriptions consider 
non-local hydrodynamic effects, albeit typically assum¬ 
ing that the filament is rigid such that elastic forces 
are ignored [UllS]. The few studies that have coupled 
long-range hydrodynamics with elasticity either assume 
small deflections m or approximate the filament as a 
network of springs [Til dS], thereby oversimplifying the 
mechanics of the problem. One exception is the study of 
buckling of a straight elastic filament loaded by viscous 
stresses m- More recently, a systematic computational 
study has been performed on a discretized model based 
on Kirchhoff’s theory for elastic rods (in the form of a 
chain of connected spheres), coupled with RFT |3]. This 
study was significant in that it was the first, to the best 
of our knowledge, to report a series of buckling instabil¬ 
ities of the flagellum that arise during locomotion and 
suggested its relevance to the biological system. More¬ 
over, it addressed the important rotation-translation cou¬ 
pling. However, recent experiments [laiiT! and simula¬ 
tions m have pointed to the oversimplifying nature of 
RFT to model propulsion in a quantitatively predictive 
manner. Therefore, there is a timely need for a descrip¬ 
tion that fully couples a geometrically nonlinear elastic 
model of the filament |18] with long-range hydrodynamic 
interactions [TOKIH], along with precision experiments for 
detailed validation. 


Here, motivated by the locomotion of uniflagellated 
bacteria, we perform a combined experimental and nu¬ 
merical investigation of the dynamics of a helical elastic 
filament rotated in a viscous fluid. Our goal is to pre- 
dictively understand the underlying mechanical instabil¬ 
ities. In our precision model experiments, we reproduce 
and systematically quantify the dynamics of the filament, 
as a function of the control and physical parameters of 
the system. In parallel, we perform numerical simula¬ 
tions that model the elastic rod using the Discrete Elastic 
Rods (DER) method [5^, coupled to a viscous drag de¬ 
scribed by Lighthill’s slender body theory (LSBT) [TU] . 



FIG. 1. (a) Experimental apparatus: a helical rod (1), is 

rotated by a motor (2), inside a glycerin bath (3), that is 
enclosed by an external water tank (4) for temperature con¬ 
trol. Two orthogonal video cameras (5,6) record the rod. (b) 
Schematic diagram of the rod. (c-d) Examples of the de¬ 
formed rod at ca = 0.6 rad/s, from both experiments and sim¬ 
ulations: (cl-c2) helical and (dl-d2) buckled configurations 
(movie in m)- See text for the properties of rod and fluid. 




2 


After validating the numerics against experiments, we 
quantify the steady state configurations of the filament 
and explore the multi-dimensional phase space of the re¬ 
sulting propulsive force. Existing data on the physical 
properties of bacterial flagella is sparse m nn, given 
the experimental challenges associated with their mea¬ 
surement. As such, we seek a dimensionless description 
that encompasses the geometric parameters of natural 
flagella, with an emphasis on the propulsive force and 
onset of buckling. The phase boundary for this insta¬ 
bility is mapped out, onto which we locate a number of 
natural bacterial systems. These results motivate us to 
speculate on the potential biological relevance of the me¬ 
chanical instabilities of rotating flagella. 

In Fig. we provide a photograph of our appara¬ 
tus. As a model for flagella, we cast a series of elas¬ 
tomeric rods with vinylpolysiloxane [22j and indepen¬ 
dently varied each of the geometric parameters (axial 
length, I, or contour length, L, helix radius, R, pitch. A, 
and cross-sectional radius, Tq, or area moment of inertia, 
/ = 7rrQ/4) and material properties (the Young’s modu¬ 
lus, E of the rod was determined by analyzing the shape 
of a suspended annulus [23]). The rod was assumed to 
be incompressible (Poisson’s ratio, ly « 0.5). During fab¬ 
rication, a polyvinyl chloride tube was wrapped around 
a cylindrical object along a helical geometry, and was 
used as a mold for the rods. The density of the rod was 
adjusted by adding iron filings (Dowling Magnets) to the 
polymer, prior to casting. Once cured and demolded, the 
filament was clamped at one end, immersed in a bath 
of glycerin (20x20x30cm^), and rotated using a step¬ 
per motor [TSj. Using digital imaging, we reconstructed 
the deformed configurations of the filament and quanti¬ 
fied its dynamics. To ensure constant and reproducible 
values for the fluid viscosity, the glycerin bath was in¬ 
serted within an external water tank to accurately con¬ 
trol the temperature within ±0.5°C (Brinkmann Lauda 
RC6). By tuning the temperature, 7.6< 0[°C] < 32.4, 
we varied the viscosity of glycerin in the range 0.50 
< /x[Pa • s] <4.45 (±0.05 Pa-s). The density of glyc¬ 
erin is Pm = 1.24 g/cm^, and despite our best effort for 
density matching, our rods had a slightly higher value 
(< 5%) than glycerin, which is however included in the 
numerics. 

For our simulations, we combined DER |20|, a robust 
and efficient computational tool for the mechanics of 
rods, and LSBT [TU], a viscous force model that accounts 
for non-local hydrodynamics. Both DER and LSBT were 
independently validated against precision experiments in 
Refs. |24| and [T2|, respectively. The helical rod is de¬ 
scribed by its centerline, parameterized by the arc-length, 
s (Fig.l^,). For the fluid loading, LSBT is used to re¬ 
late uni the local velocity, u(s), and the force per unit 


length, f(s), at each point on the rod centerline: 

u(s) = ^^±/ f(s') • j(r)ds', (1) 

J|r(s',s)|>5 

where f_L(s) = f(s) • (I — t(s) (g) t(s)) is the component of 
f in the plane perpendicular to the tangent, t(s), r(s',s) 
is the position vector from s' to s, i5 = ig the natu¬ 
ral cutoff length, and J(r) = (lf| ■*" 
een tensor. Eq. Q is then discretized and cast into a 
3N sized linear system of the form U = AF, where N 
is the number of nodes of the discretized rod [T9|. At 
each time step in DER, the viscous forces, F, are evalu¬ 
ated from the velocities, U, and the matrix A that only 
depends on the geometric configuration of the rod. To 
advance in time, we apply this external force together 
with elastic forces, update the rod configuration and it¬ 
erate. Self-contact, possible only after buckling, is ne¬ 
glected throughout, although this does not compromise 
the agreement with experiments. 

We first establish a connection with existing litera¬ 
ture for a naturally straight filament rotating in a vis¬ 
cous fluid [HI and then consider naturally curved rods. 
In Fig. [^, we present experimental photographs of un¬ 
deformed (top) and deformed (bottom) configurations, 
for three representative cases of decreasing the natural 
radius of curvature of the rod, R, while fixing its con¬ 
tour length at L = 12.00 ± 0.05 cm. These three cases 
are: i) straight rod (clamped at an angle of a = 15°, 
for consistency with Ref. m), ii) moderately curved 
rod {R/L = 0.56, a = 0) and iii) highly curved rod 
{R/L = 0.29, a = 0). All other parameters for this part 
of the study were kept fixed: tq = 1.58 ± 0.02 mm, rod 
density, pr = 1.306 ± 0.002 g/cm^, E = 1255 ± 49 kPa, 
and p = 1.32 ± 0.05 Pa - s, which ensured a Reynolds 
number < 10“^. The resulting configurations (after ini¬ 
tial transients) for these three cases are found to vary 
dramatically with R/L. 

The role of natural curvature is quantified further in 
where we plot the steady state suspended height 
(vertical distance from the clamp to the bottom of the 
rod, h) normalized by the height in the non-rotating 
case, i.e. h = h/ho, as a function of the imposed an¬ 
gular velocity, lu. Excellent quantitative agreement is 
found between experiments and simulations, with no fit¬ 
ting parameters. Given that the value of the propulsion 
force at the clamp is too low to be measured experimen¬ 
tally, we extract it from the simulations at each time 
step as Ep = e 2 )ds, where f is obtained from 

Eq. Q. In Fig. we normalize the propulsive force, 
Fp = FpL"^/{El), by the characteristic bending force in 
the rod and plot Fp versus lo. Qualitative and quanti¬ 
tative differences are observed between the three cases: 
straight, moderately curved, and highly curved rods. The 
first two undergo a shape transition at w » 0.2 rad/s. 
However, the propulsion force of the straight rod is al- 
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ways positive, whereas it is negative for the moderately 
curved rod. By contrast, the highly curved rod exhibits 
a non-monotonic behavior: h first increases to reach a 
maximum value at a; « 1.0 rad/s, where buckling occurs, 
and eventually h « 1. Since the propulsion depends on 
the deformed configuration, the resulting Fp vs. w rela¬ 
tion is markedly different from the previous two cases; Fp 
first changes sign from negative to positive at w « 0.3 
rad/s, reaches a maximum at w « 2.4 rad/s and then 
changes sign again at a; « 2.6rad/s. The coupled effect 
of curvature, flexibility and fluid forces can thus produce 
nontrivial behavior in both geometry and propulsion. 

Reassured by the quantitative agreement between nu¬ 
merics and experiments, we turn to the dynamics of he¬ 
lical filaments as macroscopic analogues of bacterial flag¬ 
ella [13 Eg. These rods were rotated in the glycerin 
bath (/r = 1.6 ± 0.05 Pa • s) at angular velocities in the 
range 0 < a; [rpm] < 8. For now, we focus on a case 
with: E = 1255 ±49 kPa, pr = 1.273 ± 0.022g/cm^, 
Z = 20 ± 0.5 cm, A = 5 ± 0.5 cm, R = 1.59 ± 0.1cm, 
and rg = 1.58 ± 0.02 mm. Under these conditions, the 
Reynolds number always remains smaller than 10“^ and 
the Stokes flow assumption is appropriate throughout. 
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FIG. 2. (a) Experimental images of (al) straight rod, (a2) 
R/L = 0.56, and (a3) R/L — 0.29 in their undeformed (top 
row) and deformed state (bottom row) at ca = 3.14 rad/s. 
(b) Normalized suspended height, h, versus angular velocity, 
uj, for experiments and simulations, (c) Simulation data of 
normalized propulsive force, Fp, versus ca. 
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FIG. 3. (a) Sequence of experimental images at a; = 0.6 

rad/s. Material properties are provided in the text (movie 
in [IS]). Time series of (b) normalized height, h{t), and (c) 
normalized propulsive force, Fp(t). (d) Normalized height, h, 
in steady state {t > 360 s) versus ca. (e) Normalized propulsive 
force, Fp, in steady state versus ca, with the shaded region 
representing the standard deviation. 


In Fig. [^, we present a sequence of experimental pho¬ 
tographs for our representative helical rod rotated at 
ca = 0.6 rad/s, starting from rest. The corresponding 
time series for the normalized suspended height, h, is 
plotted in Fig. for both experiments (solid line) and 
simulations (dashed line), with good agreement between 
the two. Any mismatch arises primarily from self-contact 
that is neglected in the simulations. Initially {t < 100 s), 
h ~ I but the conhguration eventually becomes increas¬ 
ingly distorted due to the appearance of regions of chi¬ 
ral inversion, even if the axis of the helix remains ver¬ 
tical. At later times, the rod bundles and the sus¬ 
pended height reaches an approximate steady state, with 
h ^ 0.6. The time series of the normalized propulsive 
force, Fp = FpP/(El), calculated from the simulations. 


is plotted in Fig. Concurrently with the drop in h at 
t « 150 s, Fp becomes increasingly unsteady, which we 
will show arises through a buckling instability. 

In Fig. we plot the late time average of h over 740s 
past the initial transients (t > 360 s) versus co. We find 
that h ~ 1 up to Wb = 0.51 rad/s, after which it sharply 
drops. Hereafter, we shall refer to uJi, as the critical buck¬ 
ling velocity, above which fluid loading arising due to the 
rotation causes the helical filament to buckle. The corre- 
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FIG. 4. Buckling velocity, oif,, versus (a) characteristic angular velocity, (b) normalized pitch, \/l, and (c) normalized 

helix radius, R/X. In (a), one of the four parameters {E,ro, ij.,1} was varied keeping the others fixed. In (b) and (c), the helix 
pitch and radius, respectively, were varied while fixing all the other parameters, (d) Dependence of propulsive force (color 
bar) on both the normalized pitch, X/l, and u}{R/X)’^, with m = 1.96 ± 0.05, obtained at R/X — 0.2 (see [TS]). The symbols 
correspond to (T) Caulobacter crescentus (Wild) (o), (2) Rhizobium lupini (Curly) (A), @ Salmonella (Wild) (l>), @ Rhizobium 
lupini (Semicoiled) (v); ® Escherichia coli (0), ® Vibrio alginolyticus (<l) |12l 117] . Filled (and open) symbols correspond to 
the lower (and upper) bound estimates El = 10 “^® Nm^ (and El — 10 “^^ Nm^). The errorbars represent the range in angular 
velocity, cu/(27r) = 500 ± 200 Hz. 


spending Fp is plotted in Fig.j^ as a function of w and we 
find that it increases monotonically up to Wb. Note that 
a rigid helix would yield a linear dependence between Fp 
and oj [Hj (dashed line in Fig.[^). For uj < Wf,, flexibility 
of the helical filament leads to a sub-linear Fp, when com¬ 
pared to the rigid case. For w > ojb, the average value of 
Fp drops sharply, albeit with significant fluctuations (the 
shaded region in Fig. correspond to the standard devi¬ 
ation of the averaged force value). Similar to Fp, we can 
also measure the input torque, Tp, necessary to sustain 
rotation, and this allows us to quantify the propulsive 
efficiency, which we define as 77 = Fpl/Tp. The efficiency 
remains almost constant as a function of angular velocity 
for UJ < ujb, but drops to a lower value upon buckling m- 
Also, see Ref. [3] for a characterization of propulsion as 
a function of input torque. 

We proceed to rationalize the dependence of the buck¬ 
ling velocity on the physical parameters of our system. 
The viscous force acting on the helical filament scales as 
Fy ~ flujp. Regarding the helix as an effective beam 
allows us to estimate its critical buckling load as Fy ^ 
EI/P. Instability is expected to occur when Fy « Fy, 
which yields uJb = uJbEI /where Wb = uj{X/l, R/X) 
is a dimensionless function that depends on the geome¬ 
try of the helix alone. To systematically investigate the 
dependence of Wb on the various parameters, we start 
with the values for the rod used in Figs. and [^, and 
plot oJb versus ErQ/{pLl'^), for a given geometry {X/l = 4 
and R/X = 0.32). In both experiments (filled symbols) 
and simulations (open symbols), each one of the four pa¬ 
rameters {E, tq, pL, 1} is independently varied, while fixing 
the other three. We find that the data collapse onto a 
straight line, thereby supporting our scaling analysis. To 
characterize the effect of the geometry {X/l, R/X}, we in¬ 
dependently vary the helix radius, R, and pitch. A, while 
fixing E, p. To and 1. In Fig. we plot uib as a func¬ 
tion of the normalized pitch, X/l at fixed R = 1.59 cm. 


The dependence of uJb on the normalized radius, R/X, at 
A = 5 cm is shown in Fig. From both of these plots, 
we conclude that uJb varies nonlinearly with X/l and R/X. 
These parameters must be taken into account when map¬ 
ping the results from our model system to a regime that 
is relevant to bacterial locomotion, which is address next. 

Finally, we take advantage of the efficiency of our al¬ 
gorithm to provide a biologically relevant description of 
Wb. We use the parameters of the rod used in Fig. 
except that the fluid and flagellum are assumed to be 
density matched and the axial length is increased to 
I = 0.4 m, so that tq ^ {A,R, Z}. Supported by the 
data (see [E]), we approximate Wb = M{R/X) N{X/l) 
by M{R/X) = (i?/A)“™ with m = 1.96 ± 0.05 (see 
m for details). Using this result, in Fig. [^, we con¬ 
struct a phase diagram for the propulsive force, ver¬ 
sus both a) • (i?/A)™ and X/l, where uj = ujpl^/{EI) 
is the normalized angular velocity. In Fig. |^, we also 
superpose the parameter values corresponding to bac¬ 
terial flagella of specific organisms (see caption), which 
are estimated by taking the characteristic orders of mag¬ 
nitude values for p = 10“^Pa.s, El = 10“^^ Nm^ [TT] 
(this estimate ranges from 10“^^ [26] to 10“^^ [2l])) and 
w = 102 - 103 Hz lUdllET]. Moreover, the geometric 
parameters {A/Z,i?/A} for some common bacteria (see 
caption of Fig.[^ are taken from Refs. [mini. The data 
suggests that natural flagella rotate at a rate approxi¬ 
mately within one order of magnitude of uib, where we 
have taken into account the estimated range of uj (er¬ 
rorbars) and the known uncertainty in El (rectangles). 
Note that for simplicity and generality, we ignored the 
role of the cell body, and focused on a single helical fila¬ 
ment, even though a number of the bacteria considered 
here are multi-flagellated. 

Our results raise the hypothesis that the flexibility 
of flagella imposes an upper bound on propulsive force 
through uJb, above which buckling occurs. Moreover, in 
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addition to the localized bucking that can happen at the 
hook of the flagellum |4], the reconfigurations that arise 
in the post-buckling regime of the flagellum suggest the 
possibility of a novel functional turning mechanism. This 
remains an open question, however, given the current un¬ 
certainty on the known properties of flagella. As such, 
our investigation calls for additional experimental work 
to more precisely measure the properties of natural bac¬ 
terial flagella, and more accurately image their dynamics. 

We thank Roman Stocker for enlightening discussions 
and we are grateful to the National Science Foundation 
(CMMI-1129894) for financial support. 
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